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ABSTRACT 

Motivation: Metabolic engineering aims at modulating the 
capabilities of metabolic networks by changing the activity of 
biochemical reactions. The existing constraint-based approaches 
for metabolic engineering have proven useful, but are limited only to 
reactions catalogued in various pathway databases. 
Results: We consider the alternative of designing synthetic 
strategies which can be used not only to characterize the maximum 
theoretically possible product yield but also to engineer networks 
with optimal conversion capability by using a suitable biochemically 
feasible reaction called 'stoichiometric capacitance'. In addition, we 
provide a theoretical solution for decomposing a given stoichiometric 
capacitance over a set of known enzymatic reactions. We determine 
the stoichiometric capacitance for genome-scale metabolic networks 
of 10 organisms from different kingdoms of life and examine 
its implications for the alterations in flux variability patterns. Our 
empirical findings suggest that the theoretical capacity of metabolic 
networks comes at a cost of dramatic system's changes. 
Contact: larhlimi@mpimp-golm.mpg.de, or nikoloski@mpimp-golm. 
mpg.de 

Supplementary Information: Supplementary tables are available at 
Bioinformatics online. 



1 INTRODUCTION 

Metabolic reactions play a fundamental role in sustaining cell growth 
through the import of nutrients from the environment and their 
conversion into molecules rendering a viable organism. Metabolic 
reactions do not operate in isolation; they are in fact related through 
shared meta bolites and form large-scale metabolic networks (Klan it 
and Stelling. l2003l : [Schilling et ali l2000l : ISchuster et al[ boOOh . 
Altering the structure of metabolic networks provides the means 
for modulating their capabilities. To this end, metabolic engineering 
aims at re-programming cellular networks by controlling the activity 
of biochemical reactions in order to obtain a desirable output 
and has a lready found nunierous application in biotechnology and 
medicine (iFeist et a/ll2009l : lLee et a/.Ll201ll : lMuller and Hausmannl 
l201lh . From the early attempts of engineering bacterial strains, it 
had become apparent that the biosynthetic capabilities of organisms 
are limited not only by the input in a form of nutrients but also by 
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the reaction fluxes obeying the basic physico-chemical constraints. 
Therefore, constraint-based approaches have gained considerable 
attent ion in the metabolic engineering community (jPrice et 'ali 
l2003h . 

Based on the stoichiometry of the considered reactions, flux 
balance analysis (FBA) has provided the basic framework for 
investi gating reaction fluxes in a metabolic network (E dwards and 
Palsso n.l20QQl:lKauffman et a/.L |2004 iLee et a/.Ll200^: Varma and 
Palsson J 19941) . Extensions of this framework have subsequently 
allowed systematic investigations of network modifications directed 
at enhanci ng the production o f selected targets. For instance , 
QptKnoc k (iBurgard et a/.[l2003h . OptSt rain (IPharkva et al[\2004) , 
QptR eg jPharkva and MaranasL l2006h. QptForce (R a ngana than 
et al, |201Q|) an d CASOP jHaedicke and Klami l20ld) as 
well as EMILiO (lYang et g/.L 12'oilh have facilitated prediction 
and in silico verification of network modifications, including 
insertion and deletion of reactions from other species or under- 
and over-expression of gene products. All of these alternatives 
for network modifications rely on the assumption that gene 
expression is proportional to reaction flux. The findings of 
these optimization-based approaches suggest alternative pathways 
and/or flux redirections necessary to achieve the objective of 
interest. We note that the enumerated approaches, although very 
useful in designing strategies for metabolic engineering, are 
limited to th e react i ons already catalogued in v a rious databases 
(iKarp et ali I2OOOI : iKanehisa and Gotol |200Q|: iKumar et oZl 
l2012h . The performance of the suggested strategies is in turn 
determined by the maximum achievable yield of a desired 
target which is calculated with the help of various optimization 
techniques. 

Herein, we argue that quantifying the performance of metabolic 
engineering strategies in a given metabolic network should rely on 
the maximum 'theoretically' possible product yield in the network, 
given by the optimal inter-conversion of metabolites by using any 
biochemically feasible reaction. The notion of theoretical capacity 
refers to the estimated maximum production that can be delivered 
during a given period, provided a set of inputs. To emphasize the 
analogy, in the case of metabolic networks, the inputs are represented 
by the nutrients together with the set of metabolic intermediates and 
reactions, while the maximum production to be estimated is given 
by the suitably chosen objective function. Clearly, the inputs are 
given by the existing biochemical knowledge structured in a form 
of cellular networks together with the corresponding environmental 
constraints. 
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Estimation of the theoretical capacity for large-scale metabolic 
networks is all the more pressing issue due to the incompleteness 
of biological knowledge. The gaps in the current genome-scale 
metabolic networks arising as a result of the unchara c terized 
function of 20-60% of gene products (iHanson et al 1 l201Qh 
can have p rofound effect on the outcome of optimization- based 
approaches (iFeist et a7ll2009l : [Marashi and Bockmaviil201ll) . This 
provides the first justification for using the theoretical capacity as 
a golden standard. Therefore, the maximum theoretically possible 
yield lends itself as a standard analogous to the 'theoretical 
capacity' used in finance, transportation/logistics and batteries 
research (iTodorov et all 1 19991 : [Transportation Research BoardL 
[T997h . 

In addition, given the large fraction of unknown chemical 
reactions in biological systems and the immense space of 
macromolecules potentially catalyzing chemical reactions jPobsonl 
l2004l) . the most promising targets for metabolic engineering may 
not be found among the already characterized enzymes. Instead, 
synthetic alternatives obtained through systematic screening for 
chemically feasible, but so far uncharacterized, reactions may 
offer novel promising options for metabolic engineering ( B asler 
et a/.. l2012h . This becomes particularly interesting in cases where 
the synthetic reactions can be readily decomposed into pathways 
of equivalent net conversion capability, involving only known 
enzymes. 

The article is organized as follows: we first provide a brief 
overview of optimization-based approaches for metabolic network 
analysis. We then give a detailed and formal description of our novel 
concept, called stoichiometric capacitance. Our approach is designed 
to determine the theoretical capacity of a given metabolic network 
and relies on determining a feasible chemical reaction whose 
addition to the network results in maximum product yield. Therefore, 
it can be used in designing synthetic metabolic engineering strategies 
aimed at improving the yield of any product of interest. We illustrate 
the proposed concept on the Krebs cycle with the idea of verifying 
if glu cose can be synthesized from fatty acids dde Figueiredo et all 
l2009h . In addition, we apply the proposed method to improve 
biomass yield in genome-scale metabolic networks of 10 organisms 
from different kingdoms of life. Finally, we examine the implications 
of stoichiometric capacitance for the alterations in flux variability 
patterns. 

2 METHODS 

2.1 Optimization-based approaches 

Metabolic networks are composed of metabolites inter-converted by a 
set of biochemical reactions. Based on the chosen system boundary, the 
metabolites in a given metabolic network c an be categorized as 'internal' 
and 'external' ( iHeinrich and Schusteii Il993) . The structure of a metabolic 
network is fully described by its stoichiometric matrix S eW^^^ which 
is equivalen t to the incidence m atrix of the underlying weighted directed 
hypergraph iKlamt et~all l2QQ9h . The rows of the stoichiometric matrix 
correspond to the internal metabolites and each column represents a reaction. 
At steady state, the rate of formation of every internal metabolite equals the 
rate of its consumption. This is expressed by the flux balance equation 

Sv = 0, (1) 

where v G M" denotes the vector of fluxes, also termed 'flux distribution', for 
the reactions included in the metabolic network. Aside from the flux balance 
conditions, flux distributions are also subject to other constraints related to 



thermodynamics, environmental conditions and flux capacities, which are 
often modeled as hnear inequalities 

^min^^^^max^ (2) 

where v"^^'^ and v"^^^ stand for the lower and upper bounds of the flux 
capacities of the considered reactions, respectively. 

In addition, the law of mass conservation implies that aU internal reactions 
must be stoichiometrically 'mass balanced', i.e. each stoichiometric colunm 
corresponding to an internal reaction must lie in the kernel of the mass matrix, 
M, which contains the molecular sum formulas of the internal metabolites. 
More formally, 

MS^j = 0 for all internal reactions 7 e { 1 ,...,«} , (3) 

where Mkj is the k-th atomic species (i.e. chemical element) in the sum 
formula of the l-th metabohte and S^j stands for the j-th column in the 
stoichiometric matrix corresponding to the internal reaction j. Interestingly, 
this fundamental physical constraint is rarely considered in metabolic 
network analyses, although it provides a necessary precondition for the flux 
balance equation Q. Indeed, violation of equation Js) for any metabolite 
would allow for the removal or creation of atomic species out of nothing, 
thus rendering the flux balance equation meaningless. 

The linear constraints given in equations (fl and JJ) define a polyhedron 

P = {vgM" |5v = 0, v'"'"<v<v'"^^}, (4) 

which contains all possible steady-state flux distributions. Here, we are 
interested in those flux distributions that maximize a predefined objective 
corresponding to a metabolic network function, such as growth or product 
yield. We assume that the network objective of interest is a linear function 
r defined by a vector c g R", such that c^v > 0 for all v g P. As in FBA, the 
optimal value, r*, for a metabolic network with stoichiometric matrix, S, is 
obtained by solving the following linear programming (LP) problem: 

T*=max{c^v: Sv = 0, v'"^" < v < v"^^^}. (5) 

Clearly, the set of feasible fluxes for the problem in equation Q is 
narrowed down by the constraints in equations 0 and Namely, under 
these constraints, certain reversible reactions proceed only in one direction, 
while other reactions, termed 'blocked', are unable of carrying flux. A 
reaction carrying non-zero flux is termed 'unblocked'. Furthermore, due to 
these constraints, a non-zero flux of certa in reactions implies a non-zero flux 
through other ones iBurgard et all\20o4) . In particular, several reactions can 
be essential for performing the considered metabolic capability r. 

In ana logy to essential reactions for the growth of a micro- 
organism iEdwards and Palssonl l200d) , if t*>0, an unblocked reaction 
i is termed essential for t if there exists >0 such that c^v<a^Vi for 
all veP. Let the essential reactions for r be denoted by ^. If we have 

T 

= min { ^ : v e P} for all essential reactions ie^, then the optimal (i.e. 
maximum) capability r* is given by 

r*=min{a'v/: veP mdie^}. (6) 

Therefore, increasing the maximum value of the objective t* requires the 
relaxation of the upper bound imposed by the essential reactions on the 
objective r. Such a relaxation either allows essential reactions to carry more 
flux or reduces the set of essential reactions, i.e. certain reactions are no 
longer essential for r. 

In the following, we introduce the concept of stoichiometric capacitance, 
which inherently depends on the aforementioned relaxation, and illustrate 
its applications in genome-scale metabohc networks as weU as weU-studied 
pathways. 

2.2 Stoichiometric capacitance 

The relaxation of some of the linear inequalities (|2j may lead to an increase 
of the maximal objective r* in the original network. This can be achieved 
by placing the living system in an environment which is in favor of 
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the maximization of the objective r; for instance, the hving system can 
import more nutrients. The change in the objective value r due to the 
relaxation of environm ental conditions can be as sessed using shadow prices 
in linear programming dVarma and Palssonl[l993b . Alternativelv. a significant 
improvement of t can be obtained by determining a sparse vector r e such 
that r*(r)> r* with 

T*(r) = max{c^v: Sv + r = 0, v"^"" <v <v'^^^] . (7) 

Let /(r) = {/G {1, ...,m} I r/ 7^0}. The vector r can be chosen to simulate 
a boundary reaction exporting (respectively importing) the accumulated 
(respectively depleted) internal metabolites in /(r). More importantly, if the 
vector r lies in the null space of the mass matrix M, defined in equation Js)' 
then r represents the stoichiometry of a feasible chemical conversion of the 
internal metabolites /(r), up to multiplication by scalars. This may include 
any set of feasible chemical reactions with net conversion given by r. In the 
following, we refer to r as stoichiometric capacitance. We note that if the net 
conversion of the stoichiometric capacitance is thermodynamically feasible, 
its inclusion in the considered metabolic model, specified by equation j4), 
may increase the network objective, r, given the same amount of nutrients. 

In physics and engineering, the term capacitance is classically used to 
denote the change that must be added to a system to raise the electrical 
potential by one unit. Analogously, the concept of stoichiometric capacitance 
will denote the change in the metabolic structure given by the feasible 
chemical reaction whose inclusion in the network increases the objective 
of interest. The reference point for the increase is provided by the maximum 
value of the objective in the original metabolic network. 

To investigate the extent to which we can improve the objective r by 
adding a stoichiometric capacitance r to an investigated metabolic network, 
for each internal metabolite / e { 1 , . . . , m}, we first introduce a binary variable 
Xi such that iel(r) imphes Xi = l. Let /x denote a given maximum number 
of non-zero components of r and X be a predefined upper bound on the 
flux through the added reaction. The maximum capability T*(r) and the 
corresponding stoichiometric capacitance r can be obtained by solving the 
following mixed integer linear programming (MILP) problem 

r*(r) = max c^v 



subject to: 



Sv-\-r = 0, 








Mr = 0, 


(i) 


7r<0, 


(ii) 


-Xx <r< kx, 


(iii) 




(iv) 


rj = X,- = 0, for all (' ^ £, 


(V) 


veR", reR™, xe(0,ir. 





Condition (i) ensures that the stoichiometric capacitance, i.e. the added 
reaction, is chemically feasible. In other words, it ensures that r represents 
a mass-balanced net conversion that preserves the number of consumed and 
produced atoms. The inequality in (ii) guarantees that the stoichiometric 
capacitance is thermodynamically feasible, with T denoting the vector of 
the standard Gibbs free energy of formation of the internal metabolites. 
Constraints (iii) and (iv) ensure that r involves a limited number of 
metabolites from /(r), i.e. |/(r)|</x. Finally, the optional condition (v) 
allows modeling a stoichiometric capacitance that involves only a predefined 
subset E <^{1, ...,m} of metabolites. For instance, E may be the set of 
metabolites in a particular subcellular compartment or tissue. It is clear 
that further constraints can be included in the above MILP problem to 
determine a stoichiometric capacitance corresponding to a net conversion 
of only reactions that are more likely to occur in nature. 

Although the stoichiometric capacitance represents a chemically feasible 
reaction, it may not correspond to any of the characterized enzymatic 
reactions from investigated organisms. Therefore, it is interesting to address 
the problem of determining if the stoichiometric capacitance, determined by 
solving the proposed MILP problem, can be expressed as a combination of 
a given set 71 of known enzymatic reactions. The set 71 can be defined by 



using the reactions catalogued in several existing databases of biochemica l 
reactions <Kanehisa and GotollioOotlKarp et a/ll2000l:lKumar et a/.Lr2Q12h . 

Here, we define a linear programming problem whose solution resolves 
the raised issue. Suppose, the stoichiometry of the reactions in 71 is given 
by the matrix 

where AeW^'P, BeW''^ and CeW'''^ such that A is a stoichiometric 

matrix defined in the space of m internal metabolites and ^ ^ ^ describes the 

stoichiometry of the reactions in 71 which involve additional metabolites. 
An identified stoichiometric capacitance r can be considered as a net 
conversion of the reactions in 71 if there exist two vectors a e and ^ e 
such that 

Aa+B^ = r, 

(Xi(TA^i) < 0, for all / G { 1 , . . . ,p}, 

^i{TB^i + T'C^i) < 0, for all / € {1 , . . . , ^}, 

with T' denoting the vector of the standard Gibbs free energy of formation 
of the additional m! metabolites. 

3 RESULTS AND DISCUSSION 

We implemented the proposed method for determining the 
stoichiometr i c cap acitance within the TOMLAB environment 
(lHolmstromLll999l) . 

To illustrate our novel concept, we first consider the Krebs cycle 
whose mass-balanced reactions are depicted in Figure [T] Although 
a net synthesis of glucose from fatty acids may be feasi ble in the 
entire network of human metabolism (iKaleta fl/.Ll201lh . a recent 
theoretical analysis has demonstrated t hat such a conversion canno t 
be achieved by using the Krebs cycle dde Figueiredo et a/ll2009h . 
Here, we investigate whether this conversion becomes possible 
when a stoichiometric capacitance is added to the model of the Krebs 
cycle. To this end, we focus on finding a stoichiometric capacitance 
r with the objective of synthesizing glucose, where the set /(r) of 
involved metabolites consists only of the following six metabolites: 
acetyl coenzyme A (AcCoA), coenzyme A (Co A), isocitrate (Isocit), 
(S)-malate (Mai), succinate (Succ) and water (H2O). The rows of the 
matrices M and T corresponding to these metabolites are included 
in Supplementary Table SI. 

Solving the proposed MILP problem, with S given by the 
stoichi ometric matrix of the glycoly sis/gluconeogenesis and Krebs 
cycle jde Figueiredo et all l2009h . the matrices M and T for 
all the involved metabolites, and the synthesis of glucose as a 
metabolic function r, results in the stoichiometric capacitance r 
whose equation is 

Isocit + AcCoA + H2O ^ Mai + Succ + CoA. 

Interestingly, the stoichiometric capacitance r is the sum of the 
enzymatic reactions catalyzed by the isocitrate lyase (ICL) and the 
malate synthase (MAS) from the glyoxylate cycle: 

ICL : Isocit Succ + Gly , 

MAS : Gly + AcCo A + H2O ^ Mai + CoA, 

where 'Gly' st ands for Glyoxylate. This re sult is in agreement with 
the findings in lde Figueiredo etal\ (l2009h which state that glucose 
can be produced out of AcCoA using the Krebs cycle when the 
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Fig. 1. Metabolic model of glycolysis/gluconeogenesis and Krebs cycle in humans adopted from lde Figueiredo et al\ l2009h . It has been shown that a net 
synthesis of glucose from fatty acids cannot be achieved by using the Krebs cycle. A stoichiometric capacitance with the synthesis of glucose as a metabolic 
function r is colored in red. Adding this chemical reaction to the network allows for a net conversion of fatty acids into glucose 



Table 1. Comparison of the increases in biomass yield by using the stoichiometric capacitance (SC) 



Network 


FBA 


SC 


Reaction equation 


M. barkeri, iAF692 


0.03 


160847% 


6 H2O + 2 nac (C6H4NO2) = 2 glycogen (C6H10O5) + N2 


S. cerevisiae, iNDlSO 


0.10 


72% 


24 CO2 + 10 sbt_L (C6H14O6) = 14 glycogen (C6H10O5) + 19 O2 


M. tuberculosis, iNJ 661 


0.05 


486% 


11 H + 58 ppa (C3H5O2) = 22 glycogen (C6H10O5) + 3 ttdca (C14H27O2) 


E. coli, iJR904 


0.92 


830% 


12 CO2 + 20 H = 2 glycogen (C6H10O5) + 7 O2 


E. coli, iAF1260 


0.74 


3202% 


6 CO2 + 5 glc_D (C6H12O6) = 6 glycogen (CgHioOs) + 6 O2 


E. coli, i JO 1366 


0.98 


1714% 


12 CO2 + 9 succ (C4H4O4) = 8 aeon (C6H3O6) + 6 H2O2 


S. aureus, iSB619 


0.07 


44% 


5aizc (C9H11N3O9P) + arg_L (C6H15N4O2) = aicar (C9H13N4O8P) + citr_L (C6H13N3O3) 


H. pylori, iIT341 


0.69 


370% 


2 hom_L (C4H9NO3) = 2obut (C4H5O3) + acac (C4H5O3) + 2 NH4 


H. vulgare seeds 


2.76 


0.75% 


2 mal (C4H6O5) + 3 phpyr (C3H5O7P) = 5 CO2 + 3 e4p (C4H9O7P) 


B. subtilis 


0.40 


58% 


43 dha (C3H6O3) + 6 3dhsk (C7H7O5) + 9 pime (C7H10O4) = 39 glycogen (C6H10O5) 



FBA stands for the optimal biomass yield using flux balance analysis, SC denotes the increases in the optimal biomass yield when a stoichiometric capacitance, determined by 
solving the proposed MIL? problem, is added to the corresponding network. The last column shows the reaction equation of the corresponding stoichiometric capacitances. 



enzymes involved in the glyoxylate shunt are added. We note that 
the proposed MILP problem often does not have a unique solution, 
i.e. there may be more than one stoichiometric capacitance satisfying 
the imposed constraints. For instance, additional stoichiometric 
capacitances for the synthesis of glucose are given in Supplementary 
Table S2. 

To benchmark the proposed method, in the following we explore 
the possibility of increasing biomass production in the following 



meta bolic networks: Saccharomyces cerevisiae, iNDl 50 fPuart e 



et a/.. l2004 . Methanosarcina barkeri , iAF692 (jpeist et al 



2007), 



Mycobacterium tuberculosis, iNJ 661 dJamshidi and Palsson 
Escherichia coli, i J R904, /AF12 6 0, //0 1366 jpeist et al 
lOrth et all 1201 ll: iReed et ail 120031), Staphylococcus aureus. 



2006), 



12007 



iSB619 teecker and Palssor , l2005h . Helicobacter pylori, 
iIT341 (iThiele et a/.Ll2005h, Hordeum vulsare seed (Gra f ahrend - 
Belau et al. 120091) and Bacillus subtilis dOh et al 1 [2QQ3). 
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We constrain the stoichiometric capacitance to include up to 
four metaboUtes locaUzed in the cytoplasm of the 10 considered 
networks. The results for the increase of biomass yield due to the 
inclusion of the stoichiometric capacitance are included in Table[T] It 
is evident that the stoichiometric capacitance allows for a very small 
increase of 0.75% in the case of H. vulgare's relatively small and 
tightly constrained network. However, in the cases of M. barkeri, 
iAF 692 and E. coli, /AF1260, the maximum theoretically possible 
yield is 160847% and 3202%, respectively. 

To further investigate the effect of adding the stoichiometric 
capacitan ce to a given metabolic netw o rk, we used flux variability 
analysis (iMahadevan and Schillind. l2003h . To this end, we 
computed the flux ranges with and without adding the corresponding 
stoichiometric capacitance to the networks of E. coli, i JO 1366, 
S. cerevisiae, iNDl 50 andM. tuberculosis, iNJ661. Reactions whose 
flux is zero in all optimal FBA pathways are called 'excluded 
reactions' as they are not involved in optimizing the biomass yield. 
Other reactions, termed 'indispensable' for growth, carry non-zero 
flux in all optimal FBA pathways. Figure|2]shows the changes in the 
types of reactions due the inclusion of the stoichiometric capacitance 
in the three considered network models. 

By inspecting the panels in Figure |2] the following trends, 
obtained from the flux variability analysis, become apparent: (i) 
the number of indispensable reactions specific to the network with 
stoichiometric capacitance (in orange) is smaller compared to that of 
indispensable reactions specific to the original network (i.e. without 
the stoichiometric capacitance, in red), (ii) the number of excluded 
reactions specific to the network with stoichiometric capacitance (in 
green) is smaller compared to that of excluded reactions specific 
to the original network (in blue) and (iii) the number of reactions 
that switch from excluded in the original network to indispensable 
in the altered (in blue/yellow) is smaller than that of reactions that 
switch from indispensable in the original network to excluded in the 
network altered with the stoichiometric capacitance (in red/green); 
however, this holds valid only for E. coli and S. cerevisiae. 

In addition, we observe in Figure [2] that only few reactions 
become indispensable (numbers in orange) or excluded (numbers 
in green) and only few reactions change between these types when 
adding the stoichiometric capacitance in the three networks. On 
the other hand, some of the indispensable (respectively excluded) 
reactions in the original network become no longer indispensable 
(respectively excluded). Consequently, the flux of these reactions 
may vary dramatically and so affect the optimal biomass yield, 
which corresponds to the relaxation of the constraints imposed 
on the feasible flux distributions when adding the stoichiometric 
capacitance. However, we point out that by far most reactions remain 
indispensable or excluded (red/orange and blue/green intersections), 
which demonstrates that the stoichiometric capacitance may still 
strongly rely on the fluxes of reactions from the original network, 
while allowing for a significant increase in biomass yield (Table [T]). 
On one hand, this suggests that the identified stoichiometric 
capacitance allows for improving the performance of the original 
network without strongly redirecting the activity of the commonly 
used pathways, thus yielding a realistic scenario for metabolic 
engineering approaches. Nevertheless, since the changes happen 
across inspected types of reactions with and without stoichiometric 
capacitance, we conclude that the inclusion of the corresponding 
stoichiometric capacitance in each of the considered models may 
result in dramatic system's changes. 
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Fig. 2. Venn diagram of reaction types in the metabolic networks of 
E. coli, iJ01366, S. cerevisiae, iNDl 50 and M. tuberculosis, iNJ 661. The 
reaction types are determined by flux variability analysis using either the 
corresponding original FBA model (without stoi. cap.) or the altered FBA 
model (with stoi. cap.) which includes the corresponding stoichiometric 
capacitance given in Table \l\ Excluded reactions are those which are 
not involved in any optimal FBA pathway, whereas reactions termed 
indispensable for growth carry non-zero flux in all optimal FBA pathways 

4 CONCLUSION 

Our novel concept of stoichiometric capacitance allows not only 
for quantifying the theoretical capacity of metabolic networks with 
respect to a given objective but also suggests verifiable interventions 
(i.e. feasible chemical reactions) which provide this capacity, thus 
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facilitating the design of metabolic engineering strategies. The used 
optimization-based formulation takes into account physico-chemical 
constraints to restrict the search space in a biochemically meaningful 
manner. To this end, we used a mixed integer linear program that 
can be efficiently resolved with existing optimization software. It 
becomes apparent that the proposed approach can easily be extended 
to c onsider non-linear objective functions discussed in Schuetz 
et al. (l2007h . Further constraints can be included to determine a 
stoichiometric capacitance corresponding to a net conversion of only 
reactions that are more likely to occur in nature. Our theoretical 
framework was used in addressing a concrete question that has 
long been debated in biochemistry, and concerns the possibility 
to produce sugars from fatty acids. The identified stoichiometric 
capacitance nicely confirms existing experimental results. This 
example also illustrates how the stoichiometric capacitance, taken 
as an overall reaction, can be decomposed into single enzymatic 
reactions. We also provided a linear program that addresses the issue 
of determining whether the proposed stoichiometric capacitance can 
be expressed as an overall conversion of a given set of enzymatic 
reactions. Furthermore, the analysis of genome-scale metabolic 
networks from various organisms showed that the stoichiometric 
capacitance could dramatically increase biomass production. Taken 
together, we beheve that the presented approach will prove 
valuable not only in evaluating and classifying different metabolic 
networks according to their theoretical efficiency but also as a 
method of choice for designing synthetic metabolic engineering 
strategies. 
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